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Abstract 

Bilayer vesicles form readily from mixtures of charged and neutral surfac- 
tants. When such a mixed vesicle binds an oppositely-charged object, its 
membrane partially demixes: the adhesion zone recruits more charged surfac- 
tants from the rest of the membrane. Given an unlimited supply of adhering 
objects one might expect the vesicle to remain attractive until it was com- 
pletely covered. Contrary to this expectation, we show that a vesicle can 
instead exhibit adhesion saturation, partitioning spontaneously into an at- 
tractive zone with definite area fraction, and a repulsive zone. The latter 
zone rejects additional incoming objects because counterions on the interior 
of the vesicle migrate there, effectively reversing the membrane's charge. The 
effect is strongest at high surface charge densities, low ionic strength, and 
with thin, impermeable membranes. Adhesion saturation in such a situation 
has recently been observed experimentally [H. Aranda-Espinoza et al, Science 
285 394-397 (1999)]. 
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I. INTRODUCTION 



The self-assembly of colloidal particles offers an attractive route to the synthesis of highly 
ordered, nanostructured materials. Typically these materials have been extremely soft, being 
stabilized by entropic effects. 

For example, classical colloidal crystals are three-dimensional arrays of mutually repelling 
spheres Entropic effects maintain their crystalline order in spite of a density well below 
that of close-packing. As a result, these arrays are easily disrupted by small mechanical 
shear, dilution, etc. More recently, depletion forces have been harnessed to assemble spheres 
into crystalline arrays on the walls of their container [^j. Again the physical forces between 
the spheres are repulsive, and again the resulting arrays are extremely soft. 

Attempts to create strong ordered materials from physically attracting components have 
generally produced instead highly disordered aggregates. Recently, however, Ramos et al. re- 
ported the observation of robust two-dimensional crystallites formed from negatively- charged 
latex spheres introduced into a suspension of bilayer vesicles H||. The membranes forming 
the vesicles consist of a mixture of positively-charged and neutral surfactants. The immense 
electrostatic attraction between the negative spheres and positive membranes led to the 
crystallites' great strength; their ordered 2d character arose via the intermediary role of the 
vesicles as templates for the initial self-assembly of the spheres. 

In this paper we develop some of the physics of the crucial intermediate step just men- 
tioned, elaborating and extending the discussion in ||. This stage begins when the latex 
spheres are first introduced to the vesicle suspension, and lasts for hours to days. Initially 
the spheres adsorb avidly onto the vesicles, and indeed many vesicles become completely 
covered with spheres. However, a significant subpopulation of vesicles content themselves 
with only partial coverage: on these vesicles the adsorbed spheres form a self-limiting 'raft'. 
Once the raft forms, no further spheres attach to the vesicle anywhere, though they are 
present in excess. Instead, particles in suspension are seen to approach, then wander away 
from, the vesicle. 

The theory of colloidal surface interactions is vast (for introductions see Our goal 

is to introduce a very simple mechanism for adhesion saturation, summarized graphically in 
Fig. H] below, then present some calculations to show how it works in the parameter regime 
relevant to experiments. We will argue that our effect should be qualitatively unchanged 
after many other surface-interaction effects are included in the analysis, but much work 
remains to be done to show this in detail. Sect. |D] sketches the physics of our mechanism. 
Sect. |TTT| begins the analysis using linearized Poisson-Boltzmann theory, considering in turn a 
series of more complicated situations. The linearized theory is familiar and helps to connect 
the analysis to the physical picture, but it proves to be inadequate for the interesting range 
of parameter values. Thus in Sect. [IV] we upgrade to the full nonlinear theory, which proves 
to be quite easy in this context. Finally we consider the effects of ion correlations, neglected 
in Poisson-Boltzmann theory, in Sect. [V} A glossary of symbols appears in the Appendix. 

II. PHYSICAL PICTURE 

We first briefly review the physical picture developed in [|| and summarized in Fig. [I]. 
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Consider first two dielectric surfaces bearing fixed charge densities o± of the same mag- 
nitude but opposite sign in an electrolyte solution. When they are separated by several 
screening lengths they feel little mutual attraction, since each maintains a neutralizing cloud 
of counterions. As the surfaces approach closer, eventually their screening clouds begin to 
interpenetrate. Then negative counterions from the positive surface, and positive coun- 
terions from the negative surface, can escape to infinity without violating overall charge 
neutrality. The corresponding gain in entropy reduces the system's free energy: counterion 
release drives the surfaces into contact. 

Next consider the case of two surfaces of opposite sign and unequal magnitude; for 
instance, suppose that a + < |a_|. In this case counterion release will be incomplete; after 
exhausting all the negative counterions, some positive ones will remain, trapped by the 
requirement of charge neutrality. The osmotic pressure of the trapped ions will prevent the 
surfaces from coming into perfect contact. If one surface has variable charge density, say a + , 
then additional surface charges will be pulled into the contact region in order to improve the 
contact with the approaching negative surface || . A surface charge density can for instance 
vary because the composition of the surface is variable: for instance, the surface may be a 
mixture of charged and neutral surfactants, as in the experiments of [f|||. In this case the 
recruitment of charge to the contact region will deplete the other regions, in turn rendering 
them less attractive to additional negative dielectric objects. Fig. [I]c depicts this situation: 
surfactant rearrangement in the outer monolayer of the membrane has permitted the release 
of two more ion pairs than would otherwise (panel b) have been possible. 

The rearrangement of membrane charges is limited: the relative concentration of charged 
surfactants cannot exceed unity. The maximum charge density on the outer monolayer may 
still be less than that of the approaching dielectric, and so the final contact may still be 
imperfect, as shown in Fig. [I]c. We will assume this to be the case in the rest of this paper. 

Nevertheless, a further reduction in free energy density from Fig. |l]c is still possible, once 
we remember that the inner membrane monolayer and its counterions need not play a passive 
role. Fig. shows how the remaining trapped counterions in panel c can leave the gap, even 
if the membrane is impermeable, by following the dashed horizontal arrows in panel c. After 
this rearrangement some of the charge on the negative dielectric is neutralized by surfactants 
on the inner monolayer, whose own interior counterions migrate to the nonadhesion region. [] 
Panel c also shows a rearrangement of the surfactants on the inner monolayer, further 
depleting the charge of the noncontact zone. 

Fig. raises an intriguing question: will the migration of interior counterions ever 
overwhelm and effectively reverse the charge of the membrane as seen from outside, as 
shown in the figure? Of course, cartoons alone will not settle this question, but we can 
argue physically that such an effect may well happen as follows. First we note that the 



Even if the dielectric's charge exceeds twice the monolayer charge density, as assumed in the text 
below, additional ± ion pairs can be brought from the membrane interior, with the positive ions 
remaining in the adhesion region to help neutralize the dielectric and the negative ones migrating 
to the nonadhesion region, driving its net charge still more negative. 
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FIG. 1. ('aj Cartoon of the situation. A large vesicle of mixed neutral and positively-charged surfac- 
tants attracts a limited number of negatively-charged spheres, then saturates. The Debye screening length, 
typically about 10 nm, is much smaller than the sizes of the objects. 

(b) Disposition of counterions when an approaching negative object (shaded, above right) is still far from 
the vesicle. The vesicle interior is at the bottom of the figure. The zeros denote neutral surfactants, plus 
signs the charged surfactants. Circled ± signs denote counterions in solution. The solid vertical lines joining 
charges are fictitious elastic tethers representing intuitively the electric field lines; the requirement of charge 
neutrality translates visually into the requirement that all charges be tied in this way. 

(c) Redistribution of charges when the negative dielectric object approaches the membrane, if we artificially 
forbid any electric field inside the membrane. Four pairs of counterions have been released to infinity (upper 
left). The interior monolayer, and its counterion cloud, are unchanged from (b). Zone "n" presents a net of 
one positive charge to the vesicle exterior and so remains attractive to further incoming negative objects. 

(d) The resulting state after we relax the constraint of zero electric field inside the membrane, allowing 
the ion migrations indicated by the horizontal dashed arrows in (c). One additional counterion pair has 
been released to infinity and the adhesion gap has narrowed. The net charge of the bilayer plus interior 
counterions in zone "n" has reversed sign relative to (c), and so this zone repels additional incoming negative 
objects. [Adapted with permission from M. © 1999 American Association for the Advancement of Science.] 
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positive charge density of the noncontact zone "n" is already very small in Fig. |T]c, since the 
interior monolayer and its counterion cloud cancel, and as we will see below the electrostatic 
interaction driving the depletion of charged surfactants from the outer monolayer is very 
strong. Thus only a small migration of interior counterions will suffice to get charge reversal. 
Second, the entropic cost of creating a nonuniform charge density in the interior counterion 
cloud is quadratic in the amount of charge which migrates, since the uniform distribution 
is an equilibrium state. But the free energy gain from this redistribution is linear in the 
amount of charge migration, being dominated by the derivative d/dcr + of the attractive self- 
energy (see formula ( |3.9| ) below). Thus a finite amount of counterion migration will occur, 
and this amount may well exceed the small net charge on the nonadhesion region, effectively 
reversing it. 

The rest of this paper is devoted to a quantitative justification of the intuitive argument 
just given. Before passing on to the analysis, we should remark on another feature of Fig. ^d. 
Charge reversal requires that electric fields (represented schematically by the vertical lines 
in the figure) penetrate the interior of the membrane. Since the membrane interior is a 
low dielectric constant medium, the energetic cost of these fields can be significant, another 
term quadratic in the amount of charge migration from panel c to d. If the membrane 
is sufficiently thick, this cost will reduce the charge migration below the point of charge 
reversal, a point we will need to examine quantitatively in Sect. |LV C| below. 



III. LINEARIZED MEAN-FIELD THEORY 

In this section we begin the mathematical implementation of the ideas in Sect. |TJ. We 
begin with the linearized (Debye-Hiickel) limit of low charge density, even though ultimately 
we will argue that the experiments studied here require a full nonlinear treatment. We do 
this partly because of the simplicity of the formulae, and partly to make contact with earlier 
work. 

To fix notation and keep the article self-contained we begin by rederiving some key results 
from 11, [J. The Appendix summarizes our units and all symbols used throughout the 
paper. 



A. Basic formulae 

The electrostatic potential energy of a distribution of free charges of density p(r) is 
\ J drp(r)-0(r), where ip is the electric potential |12|]. The potential created by a single 



point charge q in an infinite, uniform, dielectric medium is ^(r) = g/47re|r|. In a more 
complicated situation, is related to p(r') by some Green function G(r, r') and obeys 
Poisson's equation, V 2 ip = —p/e. 

We first imagine a uniform charge distribution of density a on the surface {z = 0}.[] The 
halfspace z < is filled with a dielectric with no free charges, and so the electric field must 



2 The assumption of fixed charge is appropriate for surfaces with fully-ionized groups at the pH 
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everywhere vanish here. The other halfspace is a univalent salt solution in equilibrium with 
a reservoir at concentration n. The reservoir must remain neutral, but it can supply ion 
pairs at a cost in free energy given by a chemical potential fik^T. The total free energy of 
the mobile ions near the surface is then F = F ent + F es , where the entropic and electrostatic 
energies in mean-field approximation are || 

F ent = k B T J dr { n + (lnn +Vo - 1) + n_(lnn_„ - 1) - ,(n + + n_) + £(n + ~ n. + n,)] 

(3.1) 



F. 



cs 



\ J drdr , en tot (r)G(r,r')en tot (r') . (3.2) 



In the above formulae, n± are the number densities of ions, while n tot (r) = p/e = n + — 
ri- + ri{ is the total signed density, including fixed surface charges with signed density 
rif. We introduced a Lagrange multiplier £ to enforce overall neutrality. The symbol vq 
is a microscopic volume factor which will drop out of all physical results. We have fixed 
the arbitrary constant in F es by setting the electrostatic energy to zero when the mobile 
counterions form a sheet coinciding with the fixed surface charge. Thus F es is the work 
needed to pull this sheet away from the surface, and so is a positive quantity. 
In equilibrium we have = 0. Away from the plane this fixes 

n±(r)v = e ^ (r)+5) , z > . (3.3) 

Here ip = eip/ksT and we have fixed the additive constant in ip by choosing ip(oo) = 0. 
Since n + = n_ = n at infinity, we get £ = and /z = lnn^o, or 

n±(r) = ne^W . (3.4) 

Substituting then gives the free energy 

F = kftTn J dr ipsmhip — 2cosh^ + |(nf/n)^ + 2 . (3-5) 



The last term of (|3.5|) is a constant which we have added by hand to cancel a term propor- 
tional to the volume of the world. 

Eqn. (|3.5|) simplifies if the dimensionless potential ip is everywhere <C 1; in this case we 
simply get F = k^T J dr |nf0. Since the fixed charge n/ is confined to a plane, the free 
energy is a purely surface term once ip has been found. 

To find '0, we note that it satisfies the Poisson equation, a property of the Green function 
used to define it. Using the charge density en±(r) found above in ( |3.4j ) gives the Poisson- 
Boltzmann equation, 



used, such as those in the experiments of |4|||. We also implicitly assume that the surfactants used 
are insoluble in water, so that their numbers in the membrane are fixed. This assumption may 
need further scrutiny, since in the experiments one surfactant species forms micelles. 
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V 2 ^ = ^i s inhV>. (3.6) 
ek B T 

Linearizing then gives the familiar Debye-Hiickel equation: V 2 ^ = K 2 ip, where k = 



2e 2 h/ek B T. 

The objects we want to consider are much bigger than the screening length \o = 1/k 
(see Fig. |l|a). Thus our geometry is essentially planar, and we need the planar solutions 
4>{z) = Be ±KZ to the Debye-Huckel equation. The electric field is then E = — V"0, which 
indeed decays exponentially on the length scale Ad- 

For a single wall we must choose the decaying solution to (|3.6| ). We fix the constant B by 

dz ' 



imposing Gauss's law at the surface: E = — ^-z = -z. Then B = ae/ nek B T, the solution is 



— CT6 

ip(z) = — - — —e~ KZ . (linearized approximation) (3-7) 

K€k B 1 

and the free energy per unit area of the isolated, charged surface is 

/self = F/(area) = k B TaB /2e = a 2 /2ne . (linearized approximation) (3.8) 



Another well-known solution to ( |3.6| ) arises in the opposite case of very high charge 
density, where z/5 3> 1 at the surface. In this case the Poisson-Boltzmann equation has 
a solution of "Gouy-Chapman" form: ip(z) = In 2e e 2^ T (z+Agc) 2 ' Here ^ ne ^ ree parameter 
is the offset Agc ; chosen to enforce Gauss's law: Agc = 2ek B T/ea. More highly-charged 
surfaces thus have smaller Agc and so a more nearly singular potential. The pathological 
behavior of ip at large z simply reflects the end of the regime $ ^> 1 at large enough z. Note 
that the electric field E z = 2k B T '/e(z + Agc) of the Gouy-Chapman solution is independent 
of the ambient salt concentration n, as it should be: the electric forces near a highly charged 
surface depend only on the surface charge. The salt concentration determines only the extent 
of the region in which the strong-field approximation is valid. 



B. Two dielectrics 

We minimized the free energy of an isolated surface, obtaining ( |3.8j ). To extract any 
useful work from this stored free energy, we would have to remove some constraint. One 
way to do this is to bring in another semiinfinite, planar dielectric^ bearing opposite surface 



3 Nothing is really infinite. The phrase "semiinfinite planar dielectric" will mean a finite dielectric 
object whose surface curvature is much smaller than k, whose interior contains no free charges, 
and whose volume is large enough that any interior electric field would be prohibitively expensive 
in energy. As two such objects approach, the gap t between them decreases but the total volume 
occupied by solution doesn't change; this is why the constant we subtracted from ( |3.5| ) really is a 
constant. 
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charge, thus changing the solution region from a half-space to a planar slab of thickness I. 
Let us suppose that a surface with a + > approaches another surface with er_ < 0. 

Parsegian and Gingell studied this situation in the linearized approximation fTO] , arguing 
as in Sect. [□] that the surfaces attract via counterion release until all of one species of 
counterions in the gap (the "minority" species) has been exhausted. If a + ^ W-\, a residual 
cloud of the other ("majority") species remains in the gap and the system equilibrates at a 
finite gap spacing t*. Nevertheless, the final state has less free energy per unit area than it 
did originally; the difference is the adhesion strength W. 

We could compute W by again solving a boundary- value problem as in Sect, [ill A|, but 



there is a shortcut. Suppose that o + < |er_|, so that the + counterions are the "major- 
ity" species. In mechanical equilibrium the hydrostatic pressure pushing the walls together 
vanishes. The planar Poisson-Boltzmann equation is a second-order ordinary differential 
equation, and so its solutions form a two-parameter family. One integration constant is 
fixed by the Gauss-law boundary condition on the negative wall, while in equilibrium the 
other is fixed by the condition of vanishing pressure. Hence the solution ip(z) is exactly 
the same for two walls as it is for the isolated negative wall; the only difference is that in 
the former case we truncate the solution at z = while in the latter case z extends to 
infinity. The equilibrium gap spacing is then just the value of z at which Gauss's law 
for the positive wall is satisfied: — (— ^ ) = Then ( |3.7| ) gives the equilibrium spacing 
by e Ki * = |<T_/<T + | in the linearized approximation. Note that indeed the right side of this 
formula is positive and greater than unity, as it must be since > 0. 

We now recall that the linearized approximation retains only the boundary term of ( |3.5D . 

so 

/ g ap(o + ,a-) = -7^- [o--V'(O) + <T+tl>(t*) 

= ^((^-) 2 -K) 2 )- (3-9) 
Repeating these steps for the opposite case where e+ > |er_|, we find that in general (Fig. |]o) 

/gap(o+,0--) = l/self - / S elf(0"-)| • (3-10) 

Remarkably, the simple combination formula ( p.lO|) will continue to hold in the full nonlinear 
Poisson-Boltzmann treatment of Sect. [IV A| below .Q Formula ( 3.10|) is certainly reasonable: 



when cr+ = |<r_| all counterions get released, the two surfaces coincide, and this was our 
reference state of zero energy. Also, when we reverse the signs of all the charges the free 
energy should not change; ( |3.10|) has this property. 



Finally we find the adhesion energy W as || 

W = /seif(c+) + /seif(o-_) - / g ap(c+, <?-) = min{((T + ) 2 , (oJ) 2 }/ek. (linearized approximation) 

(3.11) 



4 Behrens and Borkovec have independently used this fact to simplify the study of nonlinear PB 
solutions 
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Note that W is completely independent of the majority charge density, a property noted by 
Nardi et al. In light of the physical picture in Sect. ||, we can readily interpret that fact: 
The total counterion release is limited by the smaller of the two counterion populations. 

Since W is always positive we find, as expected, that oppositely- charged dielectrics always 
attract via the counterion-release mechanism |nj . Of course this is not the behavior we were 
seeking to explain (see Sect. |). We must now proceed to generalize the above arguments, 
incorporating the relevant differences between the above system and the one studied in the 
experiments of 0. 



C. Thick membrane 

We just found that two oppositely-charged dielectrics attract, as expected. But the 
experiments we are studying involve dielectric (latex) spheres interacting not with other 
dielectrics, but with a bilayer membrane. In this subsection we begin to incorporate the new 
physics associated with this situation. We first study the interaction of a dielectric of fixed 
charge density cr_ < with a positively-charged, very thick, membrane, recapitulating some 
results of Nardi et al. [[J. 

The new physical feature of this situation is that the bilayer membranes in the experi- 
ments are fluid mixtures of positively-charged and neutral surfactants. This means that the 
charge density a + on the membrane is not a fixed number, but may vary subject to a + > 
and the overall constraint that the total membrane charge / dA a + is fixed. Let cx +iav denote 
the average charge density, so that the total membrane charge is Aa +AY . In addition we 
will suppose that the charge density cannot exceed a maximum of cr max = 2e/ao determined 
by the area per headgroup oq of the charged surfactants in each of the two monolayers 
constituting the membrane, and that \a~\ > cr max . 

Throughout this paper we will adopt a highly simplified, generic picture of membrane 
compositional changes, retaining only the entropy of mixing of the two surfactant types. 
Thus we neglect other entropic or enthalpic packing effects in the assumed membrane free 
energy f m . Moreover, at first we will for simplicity neglect the bilayer structure of the 
membrane; later on we will use formulae appropriate to a bilayer. With these simplifications 
f m takes the form 

2 - ~ r(T+ ln^ + (l-^)ln(l ° + 



fm = —k B T 
a 



(3.12) 



As discussed in Sect, [nj we wish to explore the possibility of a spontaneous partition of 
the membrane into two uniform regions, which we will call zones "a" and "n" . (Ultimately 
we hope to find that "a" is adhering while "n" is nonadhering, but for the moment these 
are arbitrary names.) The areas of the two zones are not known in advance, but they must 
add up to the total area A, so we take them to be ^A and (1 — j)A respectively. 

The two zones exchange one conserved quantity, namely membrane charge .0 Thus the 



A second conserved quantity, the total charge of the counterions, is not independent but instead 
fixed by charge neutrality. 
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system can divide into zones of charge density a+ and 0+ , subject to 

7 ai a) + (1 - 7 )4 W) = ^+,av (3.13) 

This separation will be energetically advantageous if the corresponding total free energy 
A(^f(a^) + (1 — 'y)f(cr^)) is less than Af(a +iSCV ). Here the free energy density f(c + ) is 
computed for a uniform zone with fixed membrane charge density <r + , minimizing over all 
other variables. 

The instability just described will not occur if f(cr+) is a convex function, i.e. d 2 //dcx + 2 > 
0. If d 2 //dcr + 2 is negative anywhere within the allowed region < a + < cr max , we apply 
the Maxwell construction from thermodynamics to the graph of /. This involves drawing a 
straight line tangent to the graph and spanning the region of concavity. Let the two points 
of tangency be located at and cr+ . If the average membrane composition cr +jav lies 
between these two values, then the uniform system will be unstable to partitioning into two 
zones with compositions and a+\ 



In the case of a thick membrane, we have / = / gap (o"+) + /m(c+)- Consulting ( |3. 10| ) 



and (p. 12 ), we find that the first (electrostatic) term is destabilizing, while the second 



(entropic) term is stabilizing. For future use we introduce two convenient abbreviations, one 
parameterizing the relative strengths of the two terms, the other a dimensionless measure 
of charge density: 

(3 = 2ha /K = Ka /47c£ B , a = a/a m&x . (3.14) 

With these abbreviations we obtain 
^ 2 

/ = max \(d~-) 2 — (o"+) 2 | + ln<7 + + (1 — er + ) ln(l — er+)) . (linearized approximation) 
2en l J 

Nardi et al. pointed out that this function has an inflection point, giving a region of insta- 
bility (Fig. |2]&), when f3 < 1/2. According to ( |3.14| ), this means that either the maximum 
charge density e/a^ must be large, or else the salt concentration h very small. Substitut- 
ing some typical values for the charge per headgroup a = 0.5 nm and salt concentration 
n = 1 mM= 0.0006 nm -3 gives j3 = 0.006, well into the regime of instability. 

Though we have found an instability, two remarks limit its interest. First, we have 
insisted that < a + < cr max , so of course the charge densities and on our two 
zones are both positive: both zones are adhesive, unlike the experimental phenomenon we 
are trying to explain. 

Moreover, we found no instability at all unless the charge density cr max is quite large 
(recall also that is assumed to be even greater than this). But at such large charges 
our linearized approximation breaks down! Our calculation becomes inconsistent just as it 
gets interesting. Much of this paper is dedicated to correcting this deficiency. We ask the 
reader to suspend disbelief momentarily while we implement the physical picture sketched 
in Sect, [n] in the linearized theory, where the formulae are simple. Our claim is that the 
physical picture is robust and holds beyond this inadequate mathematical framework; we 
will support this claim by improving the calculation in Sect. [TV| - 
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FIG. 2. (a) Sketch of the electrostatic part of the free energy (/ gap , formula ( 3.1 0| ) ) , for a thick, mixed, 
membrane approaching a dielectric. We show the case where \<J_ \ > cr max . 



(b) Sketch of the total free energy (/ gap +/m, including formula ( 3.12 )), illustrating the Maxwell construction. 
The features of the curve have been exaggerated for clarity. The membrane will partition into a highly 
attractive region with charge density and a somewhat-attractive region with u^. 



D. Thin, permeable membrane 

The previous subsection found that a highly-charged, thick, membrane can partition into 
a zone of strong adhesion and a second zone of weaker adhesion. In this subsection we will 
introduce another element of realism by accounting for the interior charges in the vesicle. 
To highlight the key role of the membrane as a barrier to counterions, we will first study the 
simpler case of a thin membrane permeable to ions, finding uniform attraction. This sets 
the stage for the more interesting case of an impermeable membrane in Sect. llTV^ below. 

Thus the new feature introduced in this subsection is that a membrane separates the 
world into two compartments, with electrolyte solution on each side (Fig. ^a). We continue 
to neglect the internal structure of the membrane, treating it as a single thin sheet of charge; 
in Sect. [IV C| below we will improve the analysis to include the bilayer structure and finite 
internal capacitance of real membranes. 

To organize the calculation we first note that once again there is only one independent 
conserved quantity exchanged laterally between zones on the membrane, namely a + . Let 



dz e{n + (z) — n_(z)) 



(3.15) 



be the areal density of mobile interior counterions. Note that unlike cr+, which must be 
positive and less than cr max , the interior density a in can in principle have any sign and 



magnitude. We will hold a m fixed while optimizing over the gap spacing I as in Sect. U IB 



above. Since in this subsection we are assuming a permeable membrane, we then minimize 
over <T in as well to obtain f(cr + ). We will suppress explicit mention of the dependence on 
the dielectric charge <r_, because <r_ is fixed. 

The free energy density / can be regarded as an interior term from fl3.8|) , plus a gap 
term, / gap from ( |3.1(J| ), plus the membrane free energy f m from ( |3.12| ). The interior term 
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FIG. 3. ( a) Schematic for the approach of a dielectric to a thin, permeable, membrane. Negative values 
of z in the text refer to the interior of the vesicle. 

(b) Sketch of the free energy density / as a function of <r m , holding a+ fixed. 



IS /self(- CTin) = {o 

so Q3.1GQ gives / f 



ap 



2 /2eK. The gap sees opposing charge densities of er_ and (cr + + o - ;, 
= |cr_ 2 — (er+ + 0"i n ) 2 |/2eK and the equilibrium spacing is £*(er in , a + ) 



The membrane free energy f m is independent of cr in , so minimizing over <j m gives (Fig. |3|b) 
°"in,* = |c— | — cr+, a positive value corresponding to = 0: the membrane comes into tight 
contact with the dielectric. Evaluating the free energy at this point gives 



2ek 



(7-i 



p (a+ In o-+ + (1 - ct+) ln(l - a+)) 



(3.16) 



Computing the second derivative we see that this time every term of / is separately 
convex: there is no instability. Our result is physically reasonable. As the positive mem- 
brane approaches the dielectric, the latter's negative (minority) counterions and some of 
the positive (majority) counterions get released to the exterior. Since we have assumed the 
membrane is permeable, the remaining positive counterions pass through it, where many 
more pair up with interior negative ions from the membrane and get released to the interior. 
Since no counterions need to remain in the gap, we get tight contact between membrane and 
dielectric, which in effect become a single object of reduced charge density cr + — |er_|. The 
electrostatic self-energy of this composite object is a convex function, entropy never favors 
phase separation, and so there is no instability. 



E. Thin, impermeable membrane 

Previous subsections have shown that charge mobility alone can lead to an instability, 
but not to charge reversal (Sect. |III C|) , and that introducing a coupled interior compart- 
ment alone does not even lead to instability (Sect. [II ID). Surprisingly, in this section and 
Sect. [IV B| below we will find that combining these two unpromising ingredients with the 
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hypothesis of a membrane impermeable to ions can lead to a charge-reversing instability.^ 
The key observation is that an impermeable membrane has two conserved quantities inde- 
pendently exchanged between zones: the membrane charge a + and the interior counterion 
charge density <7 in .[] In the previous subsection cr in could relax within a zone by passing 
through the membrane, and so we simply optimized / over it before applying the Maxwell 
construction. For an impermeable membrane we must instead apply the Maxwell construc- 
tion to both a+ and a m jointly. 

The geometry is the same as Fig. |3]a. It will shorten some formula? to define the total 
charge density 

Ot = 0+ + Oin . (3.17) 

The free energy density is the same as in Sect. [Ill JJ| , but this time we need a more general 
formulation than ( |3.16| ), since we are not simply evaluating at the optimal value of cx in . In 
fact, there are three physically separate cases we must distinguish: 

i) The membrane plus its trapped interior counterions may have greater charge density than 
the dielectric: a t > |er_|. 

ii) The membrane plus its trapped counterions may have lower charge density than the 
dielectric, but still be positive: < a t < |<r_|. 

Hi) The trapped counterions may overwhelm and effectively reverse the charge of the mem- 
brane: <j t < 0. This is the charge-reversal we seek. In this case the equilibrium distance 
between the membrane and a negative dielectric is infinity; the membrane actually repels 
incoming negative objects. 

The total free energy density in each of these cases (and still in the linearized approxi- 
mation) now reads 



'Vax 2 f/-_ x2 , f IK) 2 - (^-) 2 | if er t > 



2 KG 



s Ot - cx+) 2 + < y ;' a /I v' 2 



(a t ) 2 + (a_) 2 if a t < 



+/3(cr + lncr + + (1 — er + ) ln(l — cr + )) . (linearized approximation) (3.18) 

The function / defines a surface over the (er + , d~ t )-plane. If this surface is everywhere 
convex-down then there is no instability. If not, then it may be possible to bring a straight 
line up to the surface from below, touching it at two points of tangency but lower than 



We need not assume the membrane to be impermeable to water; because the bulk salt concen- 
tration is assumed the same on both sides, there will be no net osmotic flow. 

7 A third exchanged quantity, the net counterion charge density outside the membrane, is then 
fixed by charge neutrality: oq = —{o\ + <r+ + a J). Similarly the density of neutral surfactants 
in the membrane is not independent, being given by (2/ao) — (<r+/e). The numbers of individual 
counterions of each species are not, however, conserved, since neutral ± pairs can be exchanged 
with large reservoirs (the bulk solution inside and outside the vesicle) without macroscopic charge 
separation. 
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the surface at some point (o"+, av , 0t,av) lying between those tangency points. In this case a 
homogeneous system with average composition (a + av , a t av ) will be able to reduce its free 
energy by partitioning into zones whose compositions are given by the two points of tangency. 

We will assume that initially, when the membrane vesicle was formed and no dielectric 
spheres were present, the ions were in equilibrium across the membrane, so that half of 
them got trapped inside: <r t ,av = |^+,av We will also for illustration generally take the 
membrane composition to be half charged and half neutral surfactants, so that the mole 
fraction a +jav = |. Finally we assume the approaching dielectric to have greater charge 
density than the maximum possible value for the membrane. For illustration we take cr_ = 
—3/2. Summarizing, we will consider the illustrative case 

o+,av = | , <7t,av = !<7+,av , O- = -3/2 . (illustrative case) 

In general there may be many lines in the (0 + , at) -plane, all passing through the point 
(o+,av, |o+,av) and all exhibiting the instability. In this case we must examine all the lines and 
choose the one which gives the absolute minimum in free energy. To do this systematically, 
we label the lines by the real number p and write each parametrically as 

(a+,a t ) = (s, i^a+, av + f s) , 0<s<l. (3.19) 



In principle one could now plot / from ( |3.18| ) along the family of lines defined by ( |3.19|) 



find the points of tangency, optimize over p, and finally obtain the sought instability and 
the charge densities a[ a ^ and (7 t in the two zones as points of tangency, as described in 



Sect. [Ill C[ above. If one of these (conventionally <7 t ) proves to be negative, then we conclude 



that the membrane exhibits a charge-reversal instability, as was to be shown. In fact these 
steps are now rather easy to complete. But we have already remarked that the linearized 
theory is not accurate in the regime of high charge densities of interest to us. Accordingly 
we will now improve our theory by solving the full nonlinear Poisson-Boltzmann equation, 
then carry out the steps just described. 



IV. NONLINEAR POISSON BOLTZMANN THEORY 

A. Basic formulae 

We introduce the useful new variable 

C = e KZ . (4.1) 
It will also be convenient to define another nondimensional form of the charge density by 

a = an/2ne = 2a/ (3 . (4.2) 
and a nondimensional form of the free energy density by 

/ = 7^/. (4.3) 
nk B T 
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The general solution to the Poisson-Boltzmann equation can be written in terms of 
elliptic functions (see for example flT3l). Fortunately, however, we need only the zero-pressure 
solutions, corresponding to two walls which are free to adopt their equilibrium spacing £*, 
and these solutions consist of elementary functions |T^,|6|: 

^± = ±21n^. (4.4) 
At ( — > oo we then have ip± — > ±4/£ = ±4e~ KZ , precisely the weak-field solution we found 



in Sect. Ill A 



The Poisson-Boltzmann equation is second-order, and so its general solution has two 
integration constants. We have fixed one of these by restricting to the zero-pressure case. 
The other one enters ( |4.4| ) rather trivially, due to the translation invariance of the PB 
equation: in fl4.4|) we are free to shift z, or equivalently multiply ( by an arbitrary constant, 
thus obtaining a one-parameter family of zero-pressure solutions. In practice we will use 
( |4.4|) in the unshifted forms given above, but select the region d < ( < ( 2 to enforce Gauss's 
law at each of the two charged surfaces. 

For example, for an isolated surface of charge density cr + > we choose the solution ip + 
with one limit at infinity and the other at (+, which we choose by requiring 



a + k B T d - 
— = — ^ 

e e az 

or 



2+ 



C+=g-(l + Vl + ^/4) • (4-5) 



The free energy formula analogous to (|3.8|) is then obtained by substituting (|4.4| ) into Q3.5|) , 

to get 

-8 (2 - C+ A + i 
/ sc i f = i 2 — + 2cr+ In (nonlinear theory) . (4.6) 

Two surfaces of the same sign charge will repel to infinite separation, so we use this formula 
for each one separately. For a negatively-charged surface we simply replace a + by |<7_| in 

It is instructive to compare (|4.6| ) to the corresponding formula in the linearized approx- 
imation, formula ( |3.8| ) (Fig. |j). While the two formula agree at low charge density, the 
linearized formula overestimates the free energy by almost an order of magnitude at the 
high charge densities of interest to us. The nonlinear PB equation also predicts a narrower 
cloud of counterions than the linearized approximation at any given charge density. 

For two oppositely charged surfaces at their equilibrium spacing we can generalize the 
argument given in Sect. [Ill B| above, again obtaining ( |3.10|) . Again suppose first that a + < 



|er_|, and so ( + > By the same logic as in Sect. |III B| , the potential ip in the gap is just 
the same as that of an isolated surface of charge density cr_, but truncated at some finite 
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FIG. 4. The self-energy density of a charged surface / se if = Kfaeif/ksTn as a function of the dimension- 
less charge density a. Top curve: Linearized approximation, eqn. (|3.8|). Bottom curve: Poisson-Boltzmann 



solution, eqn. fll.q). / se if is positive in our conventions; see Sect. Ill A 



Let us abbreviate the local integrand in (|3.5| ) by $ = v/jsmhip — 2 cosh-0 + 2. This is the 
same for either of the two solutions ip±. We thus get 

/ sclf (a ± ) = ±^ + (c ± )+ /°°dc$ 

Zen Jc± 



and 



/ga P (^ + , cr_) = JL.( ff _$(C_) + a + ^(C+)) + [ C+ dC $ • 
2en JC_ 

But the last expression just equals / se if(cr_) — / se if((7+). Repeating for the opposite case 
a + > \cr_\, we get the desired combination formula flbMUp . 



B. Thin, impermeable membrane 



Proceeding now as in Sect. PTE] , we combine (PD, Q, (Q, Q, (|3TTq) , and 



( [3.14] ) to obtain the analog of the linearized formula ( |3.18|) : the nondimensional free energy 
density of the membrane+dielectric system at its equilibrium spacing t, as a function of 
the local membrane charge density a + and the net charge density a t of counterions trapped 
inside the membrane vesicle, is now 

?r \ - f (7, 7, \ I J I /self (^t) ~ fselffi- 

f(a + , a t ) - / self (<x t - a + ) + j + ^ ( __ ; 

+ |(a + ln6r + + (l-6r + )ln(l-a + )) . (4.7) 

Here <r t = a + + a in as before and f se u(a) is the function defined by (}4.5|), ( |4.6| ) evaluated at 
5" = 2a/ p. Thus at low charge densities the first two terms of (|4.7f ) contain factors of 1//5 2 , 
and so dominate the last (mixing-entropy) term, just as in ( |3.18| ). 
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FIG. 5. Free energy density /(<7_|_, at) f° r a thin, impermeable, membrane. For eas- 
ier visualization we have inverted the figure, rescaled, and added a linear function, plotting 
— /(cx+, tJt)/1000 + 218.0 + 6.0er + — 11.4er t instead of /. The solid curve is the locus of points where 
er t = 0; points to the left of this curve represent charge-reversed states. The heavy dot is the point 
((T + av , (it,av) = (1/2,1/4) representing the average membrane composition chosen for our illustrative cal- 
culation. The two hills in the graph imply, via the Maxwell construction, that the system's ground state 
consists of two coexisting zones. Since furthermore the hills straddle the solid curve, one of the zones is 
charge-reversed. [Adapted by permission from © 1999 American Association for the Advancement of 
Science.! 



We can now carry out the program outlined in Sect. |IIIE| for the illustrative parameter 
values n = ImM, ao = 0.5 nm 2 , j3 = 0.006, a_ = —3/2, and a +jav = 1/2 discussed earlier. 
Fig. § shows the surface defined by (|4.7| ). For clarity we have shown — / instead of /, so 
that thermodynamic stability would correspond to an inverted bowl shape. We have also 
tilted the graph by adding a convenient linear function to — /, to highlight the saddle shape. 
The linear function was selected by trial and error. Adding it does not change the points of 
tangency between the surface and a straight line. 

The graph clearly displays the instability we were seeking. Moreover, one of the two 
hills on the surface clearly lies to the left of the line of charge reversal, {a t = 0}. To make 
this qualitative observation precise, we must now evaluate (|4.7| ) along the family of lines 
specified by ( |3.19| ), perform the Maxwell construction on each line, and choose the value of 
p whose tangent line has the lowest value of / at the point (x + av . Fig. |6] shows the result of 
this analysis for the illustrative values p = 2.4 and 3.8, and the optimal value p = 2.9. 



The figure shows coexistence between a zone with a 



(a) 



0.95, and another zone with 



a 



0.25. The latter zone thus presents total charge density a t = —0.11 to the outside 
of the vesicle. Since this is negative, this zone is charge-reversed and deserves its name 
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0.2 0.4 0.6 0.8 1 



FIG. 6. Three slices through Fig. || along lines passing through the average-composition point 
((X_l_,(Xt) = (1/2, 1/4) . The total charge density seen from outside the membrane is taken to be 
<Xt = ((1 — p)o'+,av + p«7+)/2 for various values of p. Unlike Fig. ^] we have not inverted this graph. Two 
illustrative slices (gray curves) show p — 2.4 and 3.8. The black curve, with p = 2.9, gives the most ad- 
vantageous mixed state, since its tangent line intersects o+ = cr+, a v = 1/2 at the lowest value of /; hence 
p» ~ 2.9. The dashed line shows coexistence between an adhesion zone with — 0.95, covering a fraction 

- (n) 

7* = 0.36 of the vesicle, and a nonadhesive zone with a\_ = 0.25. The heavy dot shows the point of charge 
reversal, where at = 0. Since the points of tangency lie on opposite sides of this dot, the nonadhesive zone 
indeed presents net negative charge to the outside world. Again the curves have been tilted for viewing by 
plotting //lOOO + 10.9ct+ - 15.2. 

as a "nonadhesive" zone. Indeed the effect is large: <r t is —45% as great as the charge 
(T +jav /2 = 1/4 presented to the outside world when there are no adhering dielectrics spheres. 
Recalling that 70"+^ + (1 — 7)0+ = <7 +)av , we find that the adhesion zone covers 36% of the 
vesicle. These results were announced in [§J]. 

C. Finite thickness, bilayer membrane 

While the above results are encouraging, and show the mathematical possibility of a 
charge-reversal instability, our model needs considerable refinement before we can take its 
results seriously. In this subsection we begin this task by acknowledging the bilayer character 
of the membrane and its finite capacitance, both neglected up to this point. The results in 
this section were also announced in ||. 

Instead of idealizing the membrane as a thin sheet of charge density cr + , we now regard it 
as two sheets of charge density ua + and (l—u)cr+ representing the charged headgroups of the 
inner and outer surfactant layers respectively (see Fig. p]d). These two layers of charge are 
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separated by a dielectric layer of thickness t and dielectric constant e m , creating a capacitor 
of capacitance c = e m /t per area. We will estimate c using the value 0.01 pF//im 2 typical 



for artificial bilayer membranes and make the useful abbreviation 

tKe/e m = Ke/c . 



T 



(4.8) 



Then r fa 7 at salt concentration n = ImM, or more generally r « 7Wn/mM. 

The free energy formula ( |4.7|) used in Sect. IV B| needs only two simple modifications: 

1) Since the membrane still presents charge density a + — cx t to the interior solution and a t to 
the exterior (Fig. |l|), the terms involving / se if are unchanged. Now, however, when a- m + ua + 
is nonzero (or equivalently a t + (u — l)er + ^ 0), there will be a nonzero electric field in the 
membrane's interior, with a capacitive energy cost per area of (er t + (u — \)a + ) 2 /2c. 

2) The membrane now consists of two fluid monolayers of mixed charged and neutral sur- 
factants. Each monolayer has a maximal density ^cr max = e/ao, attained when the density 
of neutrals is zero. Accordingly we replace the mixing entropy f m from (|3.12|) by 



a 



C"max/2 0"max/2 



Casting everything into the nondimensional forms defined above, we see that we must 
add to the formula ( |4.7|) for / the expression 



(4.9) 



4r 2 

— (a* + (u- l)a + f+- 2ua + ln(2wa+) + (1 - 2ua+) ln(l - 2ua + ) 

+2(1 - u)a + In 2(1 - u)a + + (1 - 2(1 - u)a+) ln(l - 2(1 - w)a H 



To use our formulae we add ( |4.9| ) to ( |4.7|) and again hold fixed the two conserved quantities 
a + and a t . As before we must optimize over all other variables, in this case just u, before 
performing the Maxwell construction as in Sect. [IV BL We can simply optimize ((4.9|) over u, 



since u does not enter ([4.71) . However, this optimization is subject to the four inequalities 
which u must obey: 

< ua + < 1/2 , < (1 - u)a + < 1/2 . 

A graphical analysis similar to the one shown in Fig. |6]now gives (Fig. |7|) that the best line 
through (o + , a t ) = (1/2, 1/4) has p w 2 (see (|3.19|) ), with points of tangency at a+ = 0.247, 
and = 0.65. Proceeding as in Sect. [IV B| gives coverage 7* = 63% at equilibrium and 
a t = —0.003, or about —1.2% of the charge density presented to the outside world when no 
dielectric spheres are present. 

We can readily understand the qualitative features of these results. Since r is large, 
electric fields inside the membrane are energetically costly and the two sides of the membrane 
are nearly independent. Thus the interior ion charge density a in remains nearly uniform, 
and hence nearly equal to — cr +iav /2, and similarly the inner monolayer charge density uo + m 
cr +>av /2. Then the total charge density cr t ~ (1 — u)a + m o"+, av (see Fig. [j]<f), and ( |3.19| ) 
requires that either u = 1/2 or p w 2. The solution u = 1/2 is unphysical; the solution p m 2 
is just what we found numerically. To reverse its charge, the membrane must allow electric 
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FIG. 7. Free energy density along the optimal line p ss 2 for a finite-thickness, bilayer membrane. The 
long tick on the abscissa is at CT +jav , the illustrative value for average bilayer concentration studied in the 
text. The point of charge reversal is a+ = 0.25, which is the second point from the left where the graph 

- (n) 

intercepts the abscissa. Since the tangent point is slightly to the left of this, at a\ — 0.247, we again find 
a slight charge reversal. Again the graph has been scaled and shifted to bring out its structure: we have 
plotted the quantity / — 12750 + 5615ct + . 



fields in its interior (see Fig. |l](f); the high cost of doing this accounts for the very sharp 
left-hand dip in Fig. [7| compared to Fig. ^ and for the greatly reduced degree of charge 
reversal in the finite-thickness case. The charge-reversal effect diminishes for larger values of 
r, as we predicted in Sect. |J According to ( |4.8|) this means the effect will disappear either 
for thick membranes or at large enough ion strength n. Numerically we find the critical 
value to be about 20 mM, roughly as seen in the experiments of @,|]. 

Though the charge-reversal effect seems small, it is enough to cause the rejection of 
additional negative dielectric spheres. To estimate the magnitude of this effect, consider 
what is needed to increase 7 from its equilibrium value 7* to 7* + 5. To do this we must 



choose new values of + and <xV^ + e^ n \ subject to the condition (|3.13|), which now 



reads 



(7* + S)(a^ + e W) + (1 - 7* - S)(a^ + e W) 



We then minimize the total free energy F over and subject to this constraint, finding 
that the increase in F when we force a nonequilibrium value of 7 is 



= k^m 1 ( Q) 

K 2 



1* 1 1* 
in in 

J (n) J (a) 



A6 2 



In this expression f" a ^ denotes ^£ 



etc., and A is the total membrane area. Bringing an 



additional 1 /zm 2 of negative dielectric into contact with a vesicle of area 47r(10 /1m) 2 gives 
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5 = 0.00083. Evaluating numerically then gives AF pa 3000/cb^, a huge barrier to adhesion, 
and similarly if we pull away 1 /im 2 of adhering dielectric. 



V. EFFECTS OF ION CORRELATIONS 

The charge density near a highly- charged surface can become so great as to invalidate the 
mean-field theory we have used so far. The resulting changes in the force between two plates 
have been the object of intense study since the discovery that the total force can become 



attractive for two like-charged plates, in the presence of multivalent counterions |16|J17| . The 
situation we will need to study will be much simpler than that one, for two reasons. First, 
we are only interested in the free energy density at zero force (i.e., equilibrium). Secondly, 
our effect has arisen already at the level of mean-field theory. Since we consider the case of 
monovalent counterions only, i. e. in the regimes of small to moderate ion interactions, cor- 
relation effects will turn out to be a modest correction to the main, mean-field, contribution. 
Thus we are in the regime opposite to that recently studied in refs. |T8]JT9[ . 



The criterion for mean-field theory to be valid is roughly that the electrostatic poten- 
tial energy of two ions at the mean ionic separation be smaller than the thermal energy: 
e 2 n± 3 /Arte < k^T, or equivalently that n±£ 3 B < 1. Any isolated surface, no matter how 
highly charged, will have a distance beyond which this criterion is satisfied, and so our uni- 
versal Poisson-Boltzmann solution fl4.4|) will be valid there. We will call the boundary of this 
region z = z s . We find z s using ( |4.4j ) and ( |3.4| ), obtaining z s = 0.28 nm for our illustrative 
case of ambient salt concentration h — 1 mM. 

Thus any isolated planar surface has exactly the same potential (equation ( J4.4| )) as any 
other, for z > z s . The charge density o~ + enters only via the location z + of the surface in 
the coordinate z. Given a surface of charge density cr+, we compute z + = K~ 1 \n( + , where 
( + is given by (|4.5| ). If z + > z s = 0.28 nm, then mean-field theory is everywhere accurate 
and there is no correlated-ion cloud near the surface. In the opposite case, that part of the 
ion cloud lying within the layer z + < z < z s will have nonnegligible correlations. Since z + 
is always positive, this layer is never any thicker than a typical ion radius, and so may be 
treated as a two-dimensional classical charged gas PC| . This approach may be regarded as 
an approximation to other, more refined, calculations (e.g. refs. |^l] , |22l ). 

The effect of correlations will be to reduce the free energy density, as ions can avoid each 
other, reducing their electrostatic self-energy. To apply the results of Totsuji, originally 
derived for use in the study of electrons adsorbed onto liquid helium |2"D|] , we need to know 
the two-dimensional density m of counterions in the correlated layer, m simply equals 
o + /e minus the total density in the uncorrelated region z > z s . Again using flL4|) and 
( p.4| ) in the latter region, we find m = (o~+/e) — 0.81 nm" 2 . If this quantity is negative 
then there simply is no correlated layer and m = 0. Defining the plasma parameter as 
T = ttitl, the correlation energy density can then be represented by the interpolation 
formula E c = mk-QTY(— 1.07 + 2 2r+i 3 )' wmcn i s approximately valid over the range < 
T < 5000 |20[ . Fig. [8] shows the resulting change in the free energy density, obtained from 
the thermodynamic formula / c ( m ) = rnk-^T J r ^" m f^ r for the correlation contribution f c to 
the free energy density. 

When two oppositely-charged surfaces face each other, we have seen how the minority 
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FIG. 8. Correction to the self-energy density /self from counterion correlations. The upper curve is the 
Poisson-Boltzmann result (lower curve of Fig. ||) ; the lower curve includes the correction f c discussed in the 
text. 

surface will be stripped of its counterions; the majority surface may, however, have a corre- 
lated ion cloud, so we add f c (m) to its free energy density. Since in the situations of interest 
to us the majority surface is always the dielectric sphere, and since m depends only on the 
surface's charge (not on the presence of the other surface), we find that the correlation cor- 
rection to the free energy of the negatively-charged dielectric surface is a constant, and may 
simply be dropped. We do need to include the correlation free energy of the membrane's 
interior surface, but for a finite-thickness membrane this too is nearly a constant, since as 
we have seen the inner monolayer's charge density u<j + deviates only slightly from cr+ )av /2. 
Finally, in conditions of charge reversal the outer monolayer becomes isolated and can have 
an ion cloud of its own. Since as we have seen the degree of charge reversal is very small, the 
density m of this last ion cloud is very small and the correlation contribution is negligible. 

We have just outlined qualitatively why counterion correlations may be expected to have 
little effect on the results given in Sect. |IV C| . Indeed, the numerically calculated graph 
analogous to Fig. [7| is not appreciably different from that graph, and we do not display it 
here. 



VI. DISCUSSION 

We have proposed a theoretical explanation for the phenomenon of electrostatic adhe- 
sion saturation observed experimentally in @||. While the experimental system has not 
been systematically explored yet, our model reproduces qualitatively the surprising the phe- 
nomenon of charge reversal and several salient experimental facts 

1) Adhesion saturation occurs only with mixed bilayer vesicles, that is, at mole ratios cr +jav 
not too close to zero or unity. 

2) It occurs only under conditions of sufficiently low salt. 

3) The saturated state has a very definite number of adhering objects (7* is fixed for each 
vesicle). 
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Our analysis has omitted many familiar colloidal-force effects. Many of these are short- 
ranged (e.g. solvation forces), weak compared to electrostatic forces (e.g. undulation re- 
pulsion), or rapidly decreasing with distance (e.g. van der Waals forces). In addition we 
have neglected all finite ion-size effects. We believe that our conclusions will be robust when 
such effects are introduced, in part because the crucial physics of charge reversal involves 
the immediate neighborhood of the left-hand dip in Fig. |7|, namely, the separation between 
the charge- reversal point and the tangency point. But the distance i* between the mem- 
brane and dielectric diverges as we approach the charge-reversal point from the right, so 
this physics is controlled by the long-distance behavior of the forces. Certainly the exact 
location of the tangent point depends on the right-hand part of Fig. [7] as well, where our 
theory is not reliable. But this dependence is small due to the sharpness of the left-hand 
dip in the free energy density. Even if the right-hand side of the graph differs from what we 
computed, there should be a range of membrane compositions a + av greater than a+ but 
low enough to be in the left part of the graph, and hence yielding the sort of zone separation 
we have studied. 

We have examined only equilibrium states. It is quite possible that the experimental 
system of is not in equilibrium, i.e. that the observed coverage 7 is less than the 
equilibrium value 7* because the last one or two balls is initially repelled by a finite free- 
energy barrier. But our goal was to understand the surprising existence of any barrier, not 
to predict a specific value for 7*, which in any case depends on the membrane composition.^ 

The analysis suggests a number of experimental tests of our mechanism. A mixed vesicle 
adhering to a charged dielectric surface || may provide a more controlled geometry than 
that of 0; in this case adhesion saturation suggests the possibility of observing an adhering, 
yet flaccid, vesicle. A more ambitious test could be arranged by washing out the exterior 
solution, replacing it by another of different ion strength but the same osmolarity, while 
pinning a single vesicle for observation with a micropipette. Our formula? generalize readily 
to the case where the ionic strength h is different inside and outside the vesicle. In this way 
may be possible reversibly to turn adhesion saturation on and off. 
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"Moreover, the observed ball coverage does not directly give 7, since the degree of each ball's 
coverage is not optically observable; see ||. 
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APPENDIX: NOTATION 



Constants 

We work in SI units. Thus the potential around a point charge q in vacuum is if>(r) = 
q/4iT6or, where eo = 9 • 10~ 12 Farad/meter. We treat water as a continuum dielectric with 
e = 80eo; inside the membrane e m ~ 2eo- The Bjerrum length in water is Is = e 2 /4nekBT; 
thus 4ii£b = 8.7 nm. 



Parameters 

We take for illustration a typical ambient salt concentration of n = 1 mM=6 • 10 _4 nm~ 3 . 
Then the inverse Debye length is k = \J 2he 2 / ek-gT = ^n/mM/ (9.8 nm). The salt concen- 
tration inside the vesicle is the same, due to osmotic clamping. 

We suppose a mixture of surfactants, which for simplicity have equal area per headgroup 
a = 0.5 nm 2 . Then <7 max = 2e/ao is the maximum bilayer charge density and the parameter 
(3 = 2nao/n = 0.006 measures the relative importance of mixing-entropy and electrostatic 
effects. 

We use a typical artificial bilayer capacitance of c = 0.01pF//im 2 , which enters only in 

combination with the membrane thickness t via r = tne/e m ~ 7^/n/mM. 

For illustration we take the experimentally-controllable mole fraction of charged surfac- 
tants to be cx +iav = 1/2 and one half of the corresponding counterions to be trapped on 
the vesicle interior, so that a tjav = <x +>av /2 = 1/4. We also take the approaching charged 
dielectric objects to have charge density 50% greater than the membrane, or a_ = —3/2. 



Variables 

We generally denote nondimensionalized quantities with a bar or tilde: thus a = cr/cr max , 
while a = <jn/2ne = 2a/ (3. Also the free energy density / gives rise to / = nf /nksT = 
f / (6 ■ 10 3 kBT/ /im 2 )y / n/mM, while the electrostatic potential if) gives if> = eip/ksT. Var- 
ious contributions to / include the mixing entropy of membrane surfactants / m and the 
correlation contribution f c . 

The charge density a of a surface determines its Gouy-Chapman length Agc = 2ek B T/ecr. 
Various charge densities in the text are defined in Fig. ^ for example a t = a in + a + . m 
denotes the 2d number density of ions in the dense correlated cloud near a surface, m in 
turn determines the plasma parameter V = £B\^7vm. 

Geometrical quantities include the gap width £, the total membrane area A, the fraction 
7 of A in the adhesion zone and its equilibrium value 7*. The distance z from a surface is 
sometimes expressed using ( = e KZ . 
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